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ABSTRACT 

The  asymmetry  in  current  density  along  the  periphery  of 
cylindrical  scattering  obstacles  of  both  finite  and  infinite  length 
illuminated  by  a  plane  wave  electromagnetic  field  is  discussed. 
In  the  Appendix  the  same  problem  for  a  spherical  scatterer  is 
is  considered. 
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MISSILE  CIRCUMFERENTIAL  CURRENT 
DENSITY  FOR  PLANE  WAVE 
ELECTROMAGNETIC  FIELD  ILLUMINATION 


Introduction 

The  study  of  shadowing  of  plane  wave  electromagnetic  fields  by  missiles  is  of 
importance  because  RF  leakage  into  a  missile  through  holes  and  slots  in  the  wall  and 
the  RF  pickup  of  cables  furrowed  into  the  skin  leading  to  interior  circuitry  depends 
on  the  amplitude  and  the  variation  of  current  density  along  the  circumference  of  the 
scattering  obstacle.  Clearly  the  azimuthal  orientation  of  the  missile  with  respect 
to  the  incident  field  has  a  bearing  on  the  RF  receiving  characteristics  of  the  slots 
and  cables,  if  indeed  electromagnetic  field  shadowing  is  significant^^  £ 


One  topic  discussed  briefly  in  this  memorandum  is  a  theory  for  shadowing  by 

-2-4  C 

tubular  cylinders  of  arbitrary  dimensions  developed  by  Kao.  As  a  practical  mat¬ 
ter,  the  coupled  integral  equations  derived  by  him  cannot  readily  be  solved  for  very 
long  cylinders  because  of  the  limited  number  of  storage  locations  in  a  computer. 

To  circumvent  this,  the  author  presents  a  theory  for  shadowing  by  a  cylinder  of 
infinite  length.  In  the  entire  discussion  of  both  finite  and  infinite  tubes,  it  is  assumed 
that  the  electric  field  is  polarized  parallel  to  the  axis  of  the  cylinder.  This  is  the 
case  of  dominate  interest  in  the  field  of  RF  hazards  to  ordnance.  But  this  specific 
polarization  is  an  unnecessary  restriction.  The  angle  the  incident  field  makes  with 
the  axis  of  the  cylinder  may  be  arbitrary.  In  particular,  the  magnelia-^JtTinajr 

polarized  parallel  to  the  conductor.  ^ — — '^TvT 
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Discussion  of  the  Shadowing  Properties 
of  the  Finite  Length  Holiow  Cylinder 


Consider  a  perfectly  conducting  thin-wall  tubular  cylinder  of  complete  gener¬ 
ality  in  both  length  and  radius,  as  portrayed  by  Figure  1.  The  axis  of  the  cylinder 
coincides  with  the  z-axis  of  a  cylindrical  coordinate  system  (r,  0,  z).  The  tube 
extends  over  the  interval  -h  <  z  <  h  and  is  of  radius  a.  The  incident  electric  field 
E*  is  directed  parallel  to  the  axis  of  the  cylinder  and  arrives  at  the  angle  0  =  n 
radians.  This  orientation  of  the  electric  field,  which  is  termed  E  polarization,  is 
the  one  of  primary  interest  in  the  study  of  missile  shadowing. 
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Figure  1. 

Finite  Tubular  Cylinder  Illuminated 
by  Electric  Field  E^.  J*  (0,  z)  and 
Jg(9,  z)  are  the  coexistent  total 
current  densities 


For  E  polarization  (as  well  as  for  H  polarization,  i,  e„  H*  parallel  to  the 

t  ^  t 

z-axis),  an  axial  current  density  J  (9,  z)  and  a  transverse  current  density  J  (0,  z) 

Z  b 

are  induced.  These  currents  are  coupled.  The  superscript  t  is  used  to  denote  the 
total  current  density,  i.e.,  the  sum  of  the  currents  on  the  inside  and  outside  sur¬ 
faces  of  the  tubular  cylinder.  As  will  be  made  plain  later  for  certain  values  of 

tube  length  and  electrical  radius  k  a,  where  k  =  2n/\  ,  the  inside  currents  are 

o  o  o 

negligibly  small.  Under  these  circumstances  the  error  is  not  great  to  consider 
the  cylinder  capped  at  the  ends  by  metal  surfaces  as  is  the  situation  in  practical 
missiles. 
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Both  components  of  total  current  density  J*  (0,  z)  and  Jq(0,  z)  exist  simul¬ 
taneously  on  a  tube  of  finite  length  when  excited  by  E1 ;  only  the  former  exists  on 
a  tube  that  is  infinitely  long.  The  physical  explanation  of  this  fact  is  not  trivial. 

Along  a  tube  of  finite  length  or  any  finite  obstacle,  standing  waves  of  current 
and  charge  obtain.  The  current  and  charge  are  of  even  and  odd  symmetry,  re¬ 
spectively.  Maximum  charge  concentrations  occur  at  the  ends  of  the  tube  but  there 
is  a  distribution  of  charges  all  along  the  tube.  The  illumination  of  the  tube  is  asym¬ 
metrical,  so  that  different  densities  of  charge  exist  at  different  values  of  0.  These 
nonuniform  concentrations  exert  forces  that  cause  the  charge  to  flow  in  the  0  direc¬ 
tion  around  the  tube,  bringing  into  existence  J*  (0,  z).  The  maximum  value  of  this 

t  ^  t 

component  of  total  current  density  is  J  (0,  ±h).  The  amplitude  of  J .(0,  z)  decreases 

H  0 

as  one  recedes  from  the  ends  of  the  tube  toward  its  center.  At  that  point  J  (0,  0)  s 

t  ^ 

0  for  reasons  of  symmetry.  Clearly  J0(0,  z)  is  oppositely  directed  in  the  regions 
0  <  z  <  h  and  ~h  <  z  <  0. 

For  an  infinite  cylinder  E*  gives  rise  to  axial  progressive  or  traveling  waves 
of  current  and  charge.  There  can  be  no  axial  standing  wave  of  charge  or  current 
because  there  are  no  discontinuities  (ends).  Accordingly,  there  is  no  source  for  the 
current  density  Jq(0,  z)  and  this  component  of  current  density  is  not  excited  with  the 
indicated  polarization  of  the  incident  field.  Evidently,  the  total  current  exists 
solely  on  the  outside  surface  of  an  infinitely  long  cylinder. 

Circumferential  symmetry  of  the  current  density  cannot  exist  on  the  periphery 
of  either  a  finite  or  infinite  tubular  conductor  because  the  excitation  of  the  cylinder 
provided  by  an  incident  plain  wave  is  asymmetrical.  However,  there  would  be  no 
shadowing  effect  if  one  could  arrange  illumination  in  the  form  of  a  cylindrical  wave 
converging  on  the  obstacle. 

In  general,  one  is  interested  in  comparing  J*  ( n ,  z)  with  J*  (0,  z)  and  J*  (tt,  z) 

^  Z  Z  0 

with  Jq(9,  z).  This  provides  a  measure  of  the  shadowing  of  the  missile.  Remember 
that  the  cylinder  is  first  illuminated  in  the  region  0  -  rr;  a  shadow  region  may  occur 
at  several  values  of  9  on  the  surface  of  the  same  cylinder  (depending  on  the  value  of 
kQa)  and  not  at  0  =  0  radians.  It  is  very  important  to  note  that  shadowing  for  a  cylin¬ 
der  of  finite  length  is  a  function  of  z;  i.e.,  different  values  of  shadowing  obtain  for 
z  »  0  compared  to  z  -  h. 


7 


0 


2—6 

Kao  has  done  some  brilliant  work  in  the  field  of  scattering  from  a  thin 

wall  tubular  cylinder  of  completely  arbitrary  dimensions  (length  and  radius),  for 

E-  and  H-polarization.  Some  of  his  results  for  E-polarization  are  shown  in 

Figures  2  through  4.  Curves  for  the  shadowing  of  an  infinite  cylinder  are  also 

given.  Evidently  the  ratio  |j*(0,  z)  |  /  |  J*  (ir,  z)  is  smaller  for  infinite  cylinder 

z  z 

shadowing  as  compared  to  finite  cylinder  shadowing  for  the  range  of  values  of  k^h 
and  kQa  considered  by  Kao.  In  this  sense  infinite  cylinder  theory  bounds  the  prob¬ 
lem  for  0  =  ir,  0. 

Tables  of  data  on  shadowing  by  finite  length  tubes  appear  in  Reference  3,  and 
Reference  2  contains  a  program  in  Fortran  4  for  the  IBM  360/65  computer  permitting 
further  computations. 


From  the  Kao  formulation  of  the  problem,  it  is  possible  to  separate  J^O,  z) 
and  J  (0,  z)  into  inside  and  outside  currents.  For  the  orientation  of  the  field 
assumed  the  longest  wavelength  of  the  mode  that  can  exist  in  the  tube  corresponds 
to  the  TMQ1  mode.  Thus,  if  k^a  <  2.407,  there  is  essentially  no  current  in  the 
interior  of  the  tube,  and  J*  (0,  z)  and  J*  (0,  z)  are  outside  currents.  This  presup- 

Z  w 

poses  that  the  inequality  h  >  5a  is  satisfied. 

If  the  incident  electric  field  is  tilted  with  respect  to  the  axis  of  the  tube  or 
H-polarization  is  under  consideration,  TE  modes  may  be  excited.  For  cut-off  of 
the  TE^  mode,  kQa  <  1.843. 


On  Estimating  the  Degree  of  Nonuniformity  of  the 
Current  Density  and  the  Total  Electric  Field 
About  the  Periphery  of  an  Infinite  Cylinder 
for  Plane-Wave  Illumination 

3 

The  largest  cylinders  investigated  numerically  by  Kao  have  an  electrical 
length  kQh  of  l.5ff,  i.  e.,  h  =  0.75\q,  and  electrical  radii  of  koa  =  1.0,  2.0,  and 
3. 0.  The  following  simple  analysis,  based  on  infinite  cylinder  theory  for 
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Figure  3. 

Total  Current  Density  J*(9,  z) 
and  J^(9,  z)  as  Functions  of  9 
for  Fixed  Values  of  z.  The 
radii  of  the  cylinders  are  ka  = 
2;  the  half-lengths  ate  (a) 
kh  =  0.5n,  (b)  kh  =  l.Ow,  (c) 
kh  =  1.  5rr.  The  incident  field 
is  E -polarized.  The  kh— » 
case  is  also  displayed,  (per 
Kao) 


(c)w>  ■  1.5 


Figure  4.  Total  Current  Density  Jz(9,  z)  and  Jg{9,  z)  as  Functions  of  0  for  Fixed 
Values  of  z.  The  radii  of  the  cylinders  are  ka  =  3;  the  half-lengths  are 
(a)  kh  =  0.5 it,  (b)  kh  =  1.0 n,  (c)  kh  =  1.5^.  The  incident  wave  is  E- 
polarized.  The  kh-* »  case  is  also  displayed,  (per  Kao) 


E -polarization,  may  be  of  utilitarian  value  in  estimating  shadowing  effects  for  mis¬ 
siles  when  k  h  »  1  and  h/a  »  1. 
o 

The  total  field  in  the  z-direction  in  the  vicinity  of  an  infinite  cylindrical  scat' 

7 

terer  illuminated  at  0  =  rr  by  a  plane  wave  is 


(1) 


t  *  |j  <k  r)H®\k  a)  -  J  (k  a)H^2){k  r)]ej 

2*  (r,  e>  -  Eo  Y  rnLn  °  ■  °-  (2)n  ° 

2  °  h'2  (k  a) 

n  o 


Since  interest  centers  in  the  case  (kQr)  «  1,  i.  e. ,  in  the  low-frequency 
region  of  the  radio  spectrum,  it  is  sufficient  to  retain  only  three  terms  in  the  above 


Using  the  relations 


J  <z)  =  (-1)  J  <z) 
-n  n 


H(2\z)  =  (-l)nH^2)(z) 
-n  n 


Equation  (1)  may  be  written 


8)  ■ 


lfJoV1Ho ’V1  - 


H(2)(k  a) 
o  o 


[j.^rmfVoa)  -  J.OXV] 

.40  ~ _  -  -  —  •  -  ~ 


cos  9} 


For  small  arguments  (k^r)  ^  1  this  expression  reduces  to 


-  E  *n(-) 


e>  =  j  ff  2  °  rl v F(r* e)  * 


where 
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F(r,  0)  =  1  - 


(5) 


In  Equations  (4)  and  (5),  y  =  1.781  so  that  ih  Y  =  c  =  0.5772  ...  is  Euler's 
constant.  Use  was  made  of  the  small  argument  form  of  the  Bessel  functions: 


JQ(z)  -  1 
J1(z)  =  z/2 


(z) 


H 


(6) 


Note  that  the  function  F(r,  0)  is  the  factor  by  which  the  principal  rotationally  sym¬ 
metrical  part  of  the  resultant  field  near  the  cylinder  must  be  multiplied  in  order  to 
obtain  the  first-order  correction  for  small  departures  from  rotational  symmetry. 


Numerical  Example 

F(ff)  and  F(0)  are  to  be  computed  0.5  cm  away  from  the  skin  of 

the  missile.  Let  f  =  8.485  MHz  and  a  =  53.377  cm.  It  follows  that 

_2 

Xq  =  35.357  m,  and  r  =  53.877  cm.  Then  kQr  =  9.  574  x  10  , 

(2 In)  in  (r/a)  =  5.9634  x  lO-3,  [l  +  j(2M  in  (2/Xk  a)]  =  1  +  jl.5733, 

2  2  ~2  ° 
and  (1  -  a  /r  )  =  1.8476  x  10  .  (In  evaluating  in  (r/a)  it  was  found 

convenient  to  use  the  approximation  In  (1  +  z)  -  z.)  Substituting  these 

values  into  (5),  |  F(0)  j  =  0.5441  and  |  F(rr)  j  =  1.3781.  Hence  the 


shadowing  ratio  is 


F(0) 

F(ff) 


=  0. 6125  for  the  total  electric  field  0.  5  cm 


away  from  the  surface  of  the  cylinder. 
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The  current  density  J'  (a,  9)  on  the  cylinder  is  given  by 
z 


J*(a,  9)  =  H*(a,  9)  . 


(7) 


Now, 


VxE  =  -jw^H  , 


(8) 


so  that 


3E'(r,  9) 

J  (a,  9)  =  -r— 

Z  JCOP 


dr 


r=a 


(9) 


Applying  (9)  to  (1)  it  is  found  that 


/(a,  9) 

Z 


2E 


k  a 
o  o 


rVne 

HlZ)(k  a) 
n  o 


(10) 


In  deriving  (10)  the  Wronskian  relation 


J  (z)y'  (z)  -  Y  (z)j' (z)  =  2fnz  , 
n  n  n  n 


(11) 


was  employed.  Note  that  this  formula  is  valid  for  all  z;  in  particular  for  z  =  koa. 

Observe,  too,  that  —  J  (u)  =  J  (u)  ^  and  ~  Y  (u)  =  Y  (u)  3^ .  The  substitution 
dx  n  n  dx  dx  n  n  dx 

wu  =  was  also  made  in  obtaining  (10).  Again,  for  (k^a)2  «  1,  (10)  may  be 

written 


jt  (a,  9)  = 

o  o 


1 


H(2)(k  a) 
o  o 


-  1 


2  cos  9 


„(2>„  > 
H.  (k  a) 
1  o 


(12) 
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Using  (6)  in  (12)  yields 


■£<*. 


2E 


0)  = 


^oVl1  +  3 


G(a,  0) 


(13) 


where 


G(a,  0)  =  1 


ir  k  a 
o 


j 


(14) 


Numerical  Example 


For  the  practical  situation  previously  described  one  need  only 
determine  rrkQa  to  complete  the  problem.  This  factor  has  the  value 
idcoa  =  0.2980.  It  follows  that  |G(0)|  =  0.8441  and  ] G ( rr) ]  =  1.  3781. 


Hence 


G<0) 


G(v) 


=  0. 6125. 


On  comparing  the  total  field  and  current  density  shadowing  ratios,  it  is  ob¬ 
served  that  they  are  the  same.  However,  as  r  becomes  larger,  less  perturbation 

of  E1  by  the  presence  of  the  scattering  obstacle  should  be  expected,  i.  e. ,  E*  —  E1 . 
z  z  z 

Shadowing  probably  is  best  defined  in  terms  of  the  current  density  on  the  periphery 
of  the  missile,  because  it  is  this  current  that  flows  into  holes  and  slots  that  may 
exist  in  the  missile  skin. 


Conclusion  and  Desiderata 


Kao  has  developed  a  theory  for  determining  the  current  density  inside  and 
outside  of  a  tubular  cylinder  of  arbitrary  dimensions  for  both  E-  and  H-polarization 
and  has  written  a  computer  program  for  studying  such  scattering  obstacles  numeri¬ 
cally.  It  has  been  pointed  out  that  if  kQa  <  1.843  the  tube  is  below  cut-off  for  both 
TE  and  TM  modes,  and  the  total  current  density  is  essentially  the  current  density 
on  the  out3ide  surface  of  the  tube.  In  such  cases,  the  tube  may  be  considered  as 


capped  at  the  ends.  The  Kao  program  fails  for  very  long  tubes  because  of  computer 
storage  capacity.  When  kQh  »  1  and  h/ a  »  1,  the  method  proposed  by  the  author 
for  investigating  shadowing,  based  on  infinite  cylinder  theory,  will  yield  the  upper 
bound  for  shadowing  in  finite  length  structures  (0  =  n,  0).  Accordingly,  it  may  be  a 
satisfactory  approximation  to  base  all  missile  shadowing  studies  on  (14)  if  the  fre¬ 
quency  is  sufficiently  low.  This  simple  formula  is  amenable  to  slide-rule  evalua¬ 
tion  for  obtaining  the  current  density  on  the  surface  of  the  cylinder  at  various  azi¬ 
muth  angles  6. 

Evidently  no  information  on  the  value  of  |j*(0,  z)/j*  (n,  z)  |  with  changes  in 

z)  | .  The  latter  cur¬ 
rents  do  not  exist  on  infinite  structures  for  E- polarization  but  could  be  obtained  for 
H -polarization  with  no  z-dependence. 

2 

It  appears  that,  if  the  Kao  computer  program  is  rewritten  for  use  on  the 
computer  available,  missile  shadowing  can  be  investigated  completely.  But  since 
interest  centers  on  missiles  having  radii  small  compared  to  the  wavelength,  it  may 
be  that  the  problem  of  shadowing  of  an  incident  plane  wave  electromagnetic  field  by 
a  missile  simply  fades  away. 


z  is  available  nor  are  any  data  on  the  values  of 


z)/JqK 


appendix 


SHADOWING  BY  PERFECTLY  CONDUCTING  SPHERES 


Perhaps  the  prolate  spheroid  more  nearly  conforms  to  the  shape  of  a  missile 

than  any  other  body  of  finite  dimensions.  The  spheroid  has  the  advantage  over  the 

8 

tubular  cylinder  of  finite  length  in  that  its  ends  are  closed.  Taylor  has  obtained 
exactly  the  total  axial  current  distribution  along  a  prolate  spheroidal  antenna.  This 
work  is  of  great  importance  in  the  field  of  radio  frequency  hazards  to  ordnance 
inasmuch  as  RF  leakage  into  a  missile  containing  an  access  door  is  related  to  the 
total  surface  current. 

Unfortunately,  it  appears  a  formidable  undertaking  to  determine  the  current 
densities  on  the  surface  of  a  spheroid,  but  this  can  be  accomplished  for  the  spe¬ 
cialized  spheroid  -  -a  sphere.  The  author  now  undertakes  to  investigate  shadow¬ 
ing  by  spheres  for  plane  wave  illumination  for  the  case  k^a  <<  1.  This  inequality 
implies  a  sphere  of  small  radius,  or  long  wavelength  of  the  incident  field. 

The  current  densities  on  a  perfectly  conducting  sphere  for  a  plane  wave 

propagating  in  the  z-direction  with  electric  field  parallel  to  the  x-axis  (Figure  5) 
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Figure  5.  Geometry  of  the  Problem  for  Determining  the  Current  Densities 
J  and  J..  Generated  by  Plane-Wave  Illumination  of  a  Perfectly 
Conducting  Sphere 


In  these  expressions,  which  are  exact  for  any  value  of  k^a,  Eq  is  the  ampli¬ 
tude  of  the  incident  field.  The  objective  is  to  determine  approximately  J0  and 

for  n  =  1.  The  result  obtained  should  be  useful  when  k  a  *-<  1. 

o 

To  evaluate  (Al)  and  (A2)  for  n  an  arbitrary  positive  integer,  the  following 
relations  are  needed: 

(A3) 
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k=  n 
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hn  (Z)  =  J  z  e  2^  k!r(n  ~k  +T)(j22) 


(A7) 


k=  o 
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For  example,  since  P*  (x)  and  P*(x)  are  known,  P*  (x)  may  be  found  from 

■j  f  o  1  6 


(A5),  and  P^  (x)  from  (A6).  From  (A7)  it  is  determined  that 


h^2\z)  =  je  jZ/z 


(A9) 


and 

h(2)(z)  =  -  -(l  -i)e"jz  (A10) 

1  z\  z/ 

With  (A9)  and  (A10)  it  may  prove  convenient  to  employ 

h<2)  (z)  +  h<2),(z)  =  (2n  +  Dz^h  (z)  (All) 

n-1  n+1  n 

to  generate  the  spherical  Hankel  functions  of  order  n. 


In  the  present  instance  x  =  cos  8  and  z  =  k^a.  It  is  found  using  (A3)  to  (A8) 

that 


Also, 


(2n  +  1)  _  _3_ 
n(n  +1)  j2 


(A13) 


Since  k  a  <<  1  the  only  terms  that  count  are  the  ones  involving  (k  a)  in  the 
°  “jk0a  ° 

denominator.  Also,  one  sets  e  =1.  With  these  substitutions,  (Al)  and  (A2) 

become 
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provided  kQa  <  <:  1 . 

Note  that  Jg  =  0  when  0  =  *72  as  it  should.  Observe  also  that  is  analogous 

to  ,1^(9,  z)  for  the  cylinder.  Both  expressions  are  essentially  independent  of  k^a 

when  k  a  <<  1. 
o 


Evidently,  it  is  no  marvel  that  two  components  of  current  density  exist  along 
cylinders  of  finite  length  and  also  on  the  surface  of  spheres. 
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